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A Method in the Calculus of Variations, 

By Eainard B. Bobbins. 



Intbodtjction. 



The Calculus of Variations originated in the discussion of the problem of 
finding the curve of quickest descent, the brachistochrone. This problem was 
proposed to other mathematicians by John Bernoulli* in 1696, and modifications 
of it have been discussed ever since. There was essentially but one method 
used by various mathematicians before Lagrange in dealing with this problem, 
a method developed most fully by Euler, who applied it to the general problem 
of minimizing a definite integral. The essentials of this method consist in 

1) Considering a definite integral as the limit of a sum and a differential 
equation as the limit of an equation in differences ; 

2) Using the fact that the first derivative of a function vanishes for values 
of the variable minimizing the function. 

Although Euler put his work in geometric form, his argument may be 
stated more tersely in analytic form as follows : 

The problem is to find a function y(x) such that y = y(x) minimizes 
f%" o F(x,y,dy/dx)dx. Divide the interval x x n into n equal parts of length A 
and use the notation Xi—x -\-iA. Then consider the problem of determining 
(2/o > 2/i > ■ • ■ ■> Vn) sucn as to minimize 

«?>(2/o,2/i, ,y n )^F(x 1 ,y 1 ,p x )A^-F(x 2 ,y 2 ,p 2 )A+ + F(x n ,y n , p n )A, 

where p t = {y t — y^/A. A set of necessary conditions is that 

^-^AF y (x f , y if Pi ) — \F p {x i+1 , y i+1 , p i+1 ) —F p (x { , y if p t )\=0, i=l,2, , n. 

Then Euler reasoned that since an integral is the limit of a sum and a differen- 
tial equation is the limit of a difference equation, the function y(x) minimizing 
the integral must satisfy 

F y (x,y(x),y'(x))-^F p (x,y(x),y'(x))=0. 



* See Ostwald's Klassiker der Exahten Wissensohaften, Nr. 46. 
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Other methods have shown the conclusion to be correct. Furthermore, Knezer* 
has proved this by a modification of Euler's method. 

When Lagrange f introduced the term "variation" and the suggestive 
notation "h" to represent it and showed that d(hf) = $(df), analytic methods 
replaced the former geometric methods, with a gain in rigor and possibility of 
further extension of the theory. However, it would seem that this turn of affairs 
probably prevented the early recognition of as close connections between the 
problem of the Calculus of Variations and the ordinary problem of maxima 
and minima as actually exist. 

The first part of this paper is a development of the well known tests for 
a minimum of a function of several variables. It is introduced here because 
the methods of obtaining these tests are to be used later. The second part may 
be thought of as a comparison of an algebraic problem of minimizing a sum 
with the transcendental problem of minimizing a definite integral, the sum and 
the integral corresponding as in Euler's investigations. The theorems of Part 
II may be stated roughly as follows : 

I. If, as A decreases indefinitely, the solution (y lf . . . . , y n - x ) of the 
algebraic problem approaches uniformly to corresponding ordinates of a curve 
y=y(x), where y(x) is continuous together with its first two derivatives in the 
interval under consideration, then this curve y=y(x) satisfies the first three 
necessary conditions $ for a solution of the transcendental problem. 

II. If the absolute minimum of the algebraic problem approaches a curve 
with a continuously turning tangent, this curve gives the absolute minimum of 
the transcendental problem if such exists. 

In Part III the converse problem is considered; it is shown that if y{x) 
minimizes the transcendental problem, then for n sufficiently large a solution 
of the algebraic problem exists as close as we please to y{x). 

In Part IV the problem of the Calculus of Variations is considered follow- 
ing formal analogies with the algebraic problem. An important characteristic 
of the method used in the algebraic problem is that the variations (>? x , . . . . , Yj n -i) 
are unrestricted. To get a similar state of affairs in the transcendental problem 
of minimizing f*%F(x, y, y')dx, we let the second derivative play the role of 
independent variable. This leads to the entrance of multiple integrals and 
shows clearly certain formal analogies between the algebraic and transcendental 
problems. 

* "Euler und die Variationsrechnung," Abhandlungen sur Gesohichte der Mathematisohen Wissen- 
schaften, XXV. 

f Ostwald's Klassiker, Nr. 47. 
<i See Bolza, "Lectures on the Calculus of Variations," pp. 101, 102. 
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I. The Minimum of a Function of Several, Variables. 

Consider q> (x 1} x 2 , . . . . , x n ) a real function of the real variables (x lf . . . . , x n ) 
of class* C" in a certain neighborhood R of the origin, \x { \ < Jc, i=l, 2, . . . .,n. 
Then for points (x 1 , . . . . , x n ) within this region f R, 

I. <p(x 1} . ...,x n )~<p(0,....,0)= 2a#«(0, ....,o) 



<=1 



+ 1=2 2a:^(0, 0)+?(J), 

^ ! 4=1 j = l 

in which we use the notation l 2 =xl-{- . . . . -\-x 2 n , (I) represents a quantity infini- 
tesimal with I, and 

tyixt, ,x n )=^{x u ,x n ); ^> ij (x 1 , ,aj=-g— !-(%, ,x n ). 

To find necessary conditions that <p(x 1} . . . ., sc n ) have a relative minimum 
at (0, . . . . , 0), we assume that (0, . . . ., 0) gives a relative minimum. Then 
for (x lf . . . ., x n ) sufficiently near the origin, 

$(0, , 0)<$(a?!, , x n ) 

and the left-hand side of I is positive or zero. In particular, 

I'. <?.(^,0, ....,0)-4>(0,....,0) = ^ 1 (0,....,0) 

+ ^a3*u(O,....,O)+0f(0 1 )><). 

In order that this should not change sign with x x , we must have ^(O, . . . ., 0) 
= 0. In this way we obtain 

1) fc(0,....,0)=0, fr(0,....,0)=0, ...., <M0,....,0)=0. 
Then I becomes 

II. <p(x 1} ....,x n )—cp(Q, ...,,0)= 2 2x^.(0, ,0)+Z 2 (Z), 

and T becomes 

ii'. ^(%,o ) .... ) o)-^(o,....,o)=I^ 11 (o,....,o)+^K). 

For Xy sufficiently small the sign on the right in II' depends on the sign of 
4>u(0, . . . ., 0), if 4>u(0, . . . ., 0) =£0. Thus as a second set of necessary con- 

* </> (Xi, . .. ., as„) i» said to be of class 0(»») in a given region if <f>, together with all its partial 

derivatives of orders 1, 2, , m, is continuous in this region. 

f See Goursat, "Cours d'analyse, t. I, p. 119. 

46 
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ditions that (0, . . . ., 0) give a relative minimum for <p(x x , . . . ., x n ) we have 

2) *„(0,....,0)>0; t = l,2,....,n. 

To examine II further, we introduce v lf . . . ., v n denned by the equations 

Xx = lv lt x 2 =lv 2 , , x n =lv n . 

Then II becomes 

^(x u . . . .,x H )-*(0, . . . .,0)=±p\ 2 2 v,v#<,(0,....,0) + (Z)~|. 

For Z sufficiently small the sign on the right is the sign of the quadratic form, 
unless this vanishes. Then as a third necessary condition that (0, ....,0) 
give a relative minimum for $(%i, . . . ., x n ) we have 

n n 

3) The quadratic form 2 2 t;,-t>^(0, ... .,0) must be positive definite 
or semi-definite. 

n n 

Since 2 #? = Z 2 , we have 2 ff=l; consequently 3) can be stated thus: 

2 2 v<v,4><,(0, ....,0) 

n( Vl , . . . ., Vn ) = i=Lfci 5 >o. 

f=i 
We note that Xl^, . . . . , v n ) =Q.(cv 1 , . . . . , cv n ). For this reason the function 

n 

X2 takes all its values on the hypersphere in w-space, 2 v\—l. On this locus XI 

reaches i,ts lower limit, i. e., takes on a minimum value, which is positive or 
zero. Let 

2 2 «««#«, (0, • • • • » 0) =L, 2 v\^M. 
f=i y=i »=i 

For values (t^, . . . ., v n ) for which XX is a maximum or minimum, » vanishes 

for i=l, 2, .... , w. For such values we have 

a) M,= L^-- =0, i = l, 2, . . . ., w. 

Differentiating, ~— = 2 2 <£>i,t> ,• , ~ — =2v if and conditions a) become 



b) 



f $u«i + $»«,+ -\-^> ln v % =^v lt 



where X is the value of L and therefore of XX at (v 1 , . . . ., t>„), i. e., either a 
maximum or a minimum of XI. The necessary and sufficient condition that 
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this last set of equations have a solution {v x , . . . . , v n ) other than (0, . . . ., 0) 
is that 

($11— ^)$12 $1» 

$21 ($22— ^) $2» 



c) 



:0. 



$»1$»2 ($»» — ^) 

The roots of this equation in h are real* and condition 3) is equivalent to 
3') Every root of equation c) is positive or zero. For let Jl=A/ be a 
root of c). For a,=/l' equations b) are consistent. Multiplying equations b) 
in succession by v x , v % , . . . ., v n and adding, we have for 2,=?,', 

2 2 $,,(0, ....,0)v i v j =X%v 2 i , 

i=l j'=l i=l 

or 

2 2 <p ij (0,....,0)v i v j 

n A * 

2«f 

«=i 

Since the left-hand side is positive or zero, we have /l/>0. Conversely, if 
every root of c) is positive or zero, so also are all extremal values of 

£l(v lt , v n ). 

To recapitulate: If (0, . .. ., 0) is a point in the interior of a closed 
region B in which q>(x x , . . . ., x n ) is of class C", necessary conditions that 
<|>(a ; i, . . . ., x n ) have a relative minimum at (0, . . . ., 0) are: 

1) g|(0 0)=0, 1=1,2,' n. 

2 > sffe/ •°)> ' i=1 -°- "• 

n n 

3) The quadratic form 2 2 $,-,(0, ....,0)x i x i is positive definite or 

*=i j=i 

semi-definite; or 

3) The system b) has no solution other than (0, . . . ., 0) for X<0 ; or 
3) The equation c) in X has no negative roots. 

Conversely, if 1), 2) are satisfied and if equations b) have no solution 
other than (0, . . . ., 0) for ^<0, then $(#i, . . . ., x n ) has a relative minimum 
at the origin. For the least root of equation c) is the least value of Q. on the 

n n n 

locus 2 v*=l; if this is positive, so also is 2 2 v,-v,-$# (0, . . . ., 0), and there- 

i=l ' i=l ?=1 

fore $(scj., . . . ., x n ) — $(0, . . . ., 0) >0 for I sufficiently small and (x u ... .,.x n ) 
Ejt(0,....,0). 

* Proved by Cauchy, "Exercises de Mathematiques," IV, p. 145. 
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II. The Algebeaic Pboblem and its Limit. 

Consider a real function of three independent real variables, F(x, y, p) of 
class C" in a certain simply connected region R of the xy plane and for all 
values of p. Let (x , y ), (x n , y n ) be two interior points of R, x <x n ; let the 
equation x -{-nA—x n , n being a positive integer, determine the quantity A, 
and let Xi^x -\-iA, i=l, 2, . . . ., n — 1. Now consider the sum 

<?>(2/i> • • • -,Vn-x)^F (x lt y lt &=2i)A+F (x 2 , y 2 ,?hZ_Mi^A 

+ ....+F(x n , y n , yn ~l n - x ) A 

V\ 1 ••■•■, y n -i being such that the points (x { , y t ) lie entirely in R. This is 
obviously a function of precisely the type considered in Part I. We propose 
to consider for this particular function the conditions for a minimum derived 
in Part I, and at the same time to recapitulate the essentials of the argument. 
For n a fixed positive integer suppose that (Y 1 , . . . ., Y n _ x ), such that 
(x { , Y t ) is an interior point of R, i = l, 2, . . . ., n — 1, gives a relative minimum 
for $(«/i, • • • ., «/„_i). Then if y\ x , ....,n n -\ are an y n — 1 rea l quantities, 
a positive number 8 can be found such that when | e | < & 

1) <t>(Y x + ey, x , . . . ., Y n _ 1 + sy !n _ 1 )-^(Y 1 , ...., Y^J^O. 

Expanding by Taylor's series with a remainder (this expansion is legitimate 
with our hypotheses on F) and using the abbreviated bracket notation, 

2) x{F(x tt Y i+evii , *W*-i+* (m-y,^) ) __ F ^ f Io I^b^ a 
= '2 i {n«^(*)+ 54 =^^(i)}A 

in which >7 =»7„=0 and (e) indicates a quantity which is infinitesimal with e. 

Since the right-hand side of equation 2) must not change sign with s, for 
any real values (vj x , . . . . , *7»_i) we must have 



3) jk^F^i) +*^-^ F p (i)} A=0. 



Then the coefficient of each 57 in 3) must vanish. We can get the coefficient of 
rii in 3) very readily by inspection, but to make the analogy to the method 
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of the Calculus of Variations more evident we use the formula for partial 
summation : 

W V) 

2 &<(«r-<*<-i)— 2 «*(&<— 6«+i)— a v -ib v +a w b w+1 . 
Since >7 =>7 W =0, this identity gives us 

i(^),, (0 =S^*^£<i±il), 

and 3) becomes 

i i .,.[f,(0- ^ (<+1 i- y ' (<) -]A=o. 

Hence in this case conditions 1), Part I, are 

F y (i)-l ^^)-F P _iil = , i = l f 2,....,«-l. 
Also in 2) the coefficient of c 2 must be positive or zero; i. e., 

in which P,*=F w (t) ; Q^F^t) ; R^F pp (i). 

Let j be a positive integer less than n, and choose r ti = 0, i =/= ;/', yj.=f=0. 
(This is analogous to a method often useful in the Calculus of Variations of 
giving the minimizing curve an increment for a very short distance). Then 
the coefficient of yi) in 4) must be greater than or equal to zero, giving con- 
dition 2), Part I, for this particular problem, 

AP,.+2<?,+^±^i>0, i = l, 2, . . . ., n-1. 

Treating the quadratic form 4) in a manner closely analogous to that 
used in dealing with the second variation in the general algebraic problem, 
Part I, we have immediately, 

a. x {**.+***=2 t *Q.+ (^N _„ 

2>7?A 
1=1 

Since the expression A is homogeneous of degree zero in the >7's, it has 
the same value for the set (r a , . . . ., »7„_i) as for (p>7 2 , . . . ., p>7»-i)- Then it 
takes all its values on the locus y\\ + . . . . + rf n _ x — 1/A, a closed locus, and 
therefore it assumes its lower and its upper limits, and its lower limit is 
positive or zero. At points where extremal values are taken we have, from A, 

" d 

2 >7 2 A=— (numerator of A) — numerator of A • 2y! { A — 0, i—1, 2, . . . ., n — 1. 



374 Bobbins: A Method in the Calculus of Variations. 

» n— 1 

For such points on 2 ri\ A = 1 (or 2 >7?A = 1, since ^ n =0), 

~— (numerator of A) = 2>^A-/L 

where 2, is the extremal value of the fraction at the point in question. Per- 
forming the indicated differentiation and using the notation u { — Ui_ x — Su t , 

v, P « &Q<+i l+1 A \A_ Q 3»7«_n ^>7j+i _ o 7 _1 o ™_1 . 

or, re-arranging, 

B. B t+1 A + "Ar"A+V~A Pi+ljn+Qi+i-A &-£-", 

i — 1, 2, . . . ., w — 1. 

If {Y 1 , . . . ., Y n _x) gives a relative minimum, system B can have no solution 
{r n , . . . .; >7„_i)^0 for a,<0. It should be remembered that >? = >?„ = . 

Summarizing: If (Y lf Z 2 , . . . ., P„_i) gives a minimum for q>{y x , . . . . ,«/„_i), 
then 

r. ^ w _^(!±^=&(!i=o f -i=i,2,....,.-i. 

II'. A 2 P,+2A-^+i? i+1 +^>0, t=l,2, ....,n-l. 

IIP. The set of equations B have no solution (^, . . . ., J7»_i)=£0 for 
/l<0 such that >7 =>7 ra =0. 

Also, as in the general case, Part I, we have as sufficient conditions for a 
relative minimum, conditions P, IP and 

III". The set of equations B have no solution (r a , . . . ., >?„_i) Ej£0, such 
that y] =vj n =pO, for a,<0. 

We wish now to see what becomes of these conditions P, IP, III' as n, the 
number of divisions of our interval, increases indefinitely, assuming that the set 
(Y lf . . . ., P„_i) approach the corresponding ordinates of a curve of class C", 
y(x) ; we hope in this way to show that y(x) satisfies certain necessary con- 
ditions for minimizing the definite integral, f^F(x, y, dy/dx)dx. 

Let x' be a fixed abscissa of our interval, x <x'<x n , and for any par- 
ticular division of our interval let j be such that x i __ x <x' <#/. It is evident 
that j increases indefinitely with n, while Lt x^—x'; this will be indicated thus: 

y=co 

Lt Xi—af. 

?l=00 ( ^=00 
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Assume that, as we increase the number of points of division x 4 , the ordi- 
nates (Y lf . . . ., Y „_-,.) approach uniformly to the corresponding values of a 
function y (x) of class C" in the interval o? <#<;r n ; i. e., given e>0, N can be 
found independent of i such that whenever the number of divisions n of the 
interval {x x n ) is greater than N, 

\Y t —y(Xi) | <e for t=l, 2, , n— 1. 

Then because of the uniform continuity of y (x) in the closed interval (x x n ), 
given s > 0, N can be found independent of x' such that whenever n > N, 
\Yj — y{oc') | <s for any x' of our interval. 

Assume furthermore that Lt (Y t — Y i _ 1 )/A=y'(x i ), and that this limit is 

71 = 00 

approached uniformly. Then because of the uniform continuity of y' (x) in 
(x x n ), Lt (Yj — Y j_ x ) / A=y' (x r ) , and this limit is approached uniformly for 

71 = 00, ^ = 00 

all values of x'. 

These assumptions, together with the assumption that F y is a continuous 
function of its three arguments for (x, y) in region R and for all values of p, 
gives us that 

Lt F,\x„ Y } , Lzgj=k \=F,[af, y(x'), y'(x')], 

and that this limit is approached uniformly for all values of x' involved. 
(When convenient, we shall denote the arguments [x, y(x), y'(x)] by [%].) 
Just as with F y , it is evident that F p , F yy , F yp , F„ , F ypx , F m , F ppx , F ppy , F ppp , 
with the arguments (j), approach uniformly to the corresponding functions 
with the arguments [x'], in the interval x Q <x'<x n . 

We shall now consider equations I' as n increases indefinitely : 
I'. AF y (i) +F p (i)-F p (i+1) =0, t=l, 2, . . . ., n-1. 

We wish to show that y(x) satisfies Euler's differential equation; i. e., that 

I. F y [x]--^F p [x]=0. 

From I' we have 

5) F,{l)-F,ti) + 2F,{i)-A=0. 

i=l 

Since Lt Xj—x', and since F y {i) approaches F y [x x ] uniformly as n increases 

n=oo, /=oo 

indefinitely, 

Lt i i 1 F y (i)-A = fiF v [x]dx. 

n=oo_, ^=oo i=l 

Also 

UF p (l)=F p [Xo\; Lt F p (j)=F p [x'}. 

n~-<*> 7i=co, /=oo 
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Since the limit of a sum is the sum of the limits, it follows from 5) that 

F p [x o ]-F p [x'}+tf o F y [x]dx=:0. 
Differentiating with respect to x', we obtain equation I for x in the interval 

x < x < x n . ■ Since F y and j- F p are continuous, it follows that equation I 

holds for all values of x in the closed interval x <x<x n . 
We may write I' in the form 

• M+l ± J ±J *■ j- l 

F y (i)-F px {j+e 1 )-Yi±±=IiF vy u-{.e,) * — _ — ^- Fpp (j+d 3 )=o, 

in which (j + t ) indicates that the arguments are taken between those indicated 
by (j) and by (i + 1). We know that Lt of each function involved in this 

n=oo, 7 = 00 

last equation exists and is approached uniformly throughout the interval with 
the exception of the fraction in the last term. Then if we assume that 
Fpp[x]=t=0 in the interval %<a?<a;„, i.e., that F pp {x, y, p)=£0 along y(x), 
we have 

Lt ^ +1 -3r, = F y [x']-F px [x']-y'(x')F py [x'] 
»=.,/=- A 2 F„[a>'] ' 

From equation I we obtain 

y (x) ¥ - m . 

Hence, 

6) Lt *L±£rlL = y*(*), 

71 = 00, ^ = 00 i-i 

and furthermore this limit is approached uniformly for x <x'<x n . 
Consider condition II', which we rewrite, 

II'. A*P j +2AQ j +R j+1 + R i >0. 

It is evident that 

Lt \A*P j +2AQ j +R i+ ,+R j \ = 2F pp [x']. 

Thus we get, as the limiting form of the second condition, 

II. F pp [x] = R[x]>0, x <x<x n . 

Assuming as above that F pp [x] =£0, we have 

II". F pp [x]>0. 

We shall now use condition III' to show that if u(x) is a solution of 

4,(u)^-(Ru')-(P-Q')u=0, 
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such that u(x ) = 0, u(x)^0 (the arguments of P, Q', R being [a?]), then u(x) 
has no root between x=x and x—x n ; i. e., u(x)=f=0, in the interval x <x<x n . 

We shall make use here of the two following theorems : * 

a) Given the system 



'ffi(*(-)£)+[^(-)-'<«)]«=o, 



(**;&) - e ° M=0 ' 



u(a) =a', 

u'(a) = a, 

in which k(x) >0, #(#)>(), fc(«) is of class C", a<a;<6 and g(x), l(x) are 
continuous in the same interval, and a, a' are constants not both zero; there 
is one and only one value of the parameter X for which the solution of this 
system has a specified number (including zero) of roots between a and b and 
is such that u(b) = 0. 

b) Given the systems 

u(a) = 0, 

d 

dx l 

u(a) = 0, 

in which x lt x 2 are functions of x of class C", and G x , G 2 are continuous in the 
interval a<rc<&; furthermore, G 2 <G X and <x 2 <x 1 . Let u-l(x) be a non- 
identically vanishing solution of the first system and u 2 (x) likewise for the 
second. If u x (x) has k roots in the interval a<x<b, then u 2 {x) has at least 
k roots in this interval and the i-th root of u 2 {x) is smaller than the i-th root 
of u-i_(x) . 

Suppose that ^(u) ~0 has a non-identically vanishing solution, u(x), such 
that w(:r )=0, w(ic')=:0 for some number x' such that o: <a;'<ic n . To prove 
our theorem we shall show that this supposition leads to a contradiction; in 
fact it leads to a solution of the difference equation system B, (yj 1 , . . . . , k„_i)^0, 
such that y; z=zyj n —0, for 2, negative. 

Making use of the above theorems a), b), our present supposition gives us 
that there is one and only one value X o <0 such that the solution of the system 

r 4>(u)+*, u=o, 

7) - u(x o )=0, 

tt'(x )—y, 

* Given in this form by Professor B6cher in lectures. See also in article by B6cher, Bulletin of 
Amer. Math. 800. , Vol. IV, p. 298, I, and p. 304, VII, and Transactions of Amer. Math. 80c, Vol. I, p. 418. 

47 
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in which y is a positive constant, has x=x n as its first root to the right of 
x=x . If we replace X, hy /l — 8, £>0, the solution of system 7), u{x, Jl — h) 
does not vanish in the interval x ^x<x n and u(x n ) is positive; for X replaced 
hy \ + 8 the solution u(x, /l +^) vanishes within the interval (x x n ), and for 8 
sufficiently small, u(x n ) is negative. 

The plan is first to show that we can get one solution of the difference 
system' B, with >?o=0, as near to u(x, X — 8) as we please and therefore with 
>7 n >0, and one as near to u(x, % -\-o~) as we please and therefore with >? tt <0. 
Then by holding n fast and varying a, between /l +<5 and X — & we shall get a 
solution of B such that yi =y} n =0, for some negative "k, which contradicts condi- 
tion HI'. 

The solution of the linear differential system 

' 4>(u) +A,m=0, 
8) \ u(x )=0, 

[ u'(x )=y 
can be obtained by the Cauchy— Lipschitz method as the limit approached by 
the solution of the linear difference equation system (using our previous 
notation, 5^=^— m~i), 



9) 



BW^+B'W^+CC'W-PWI^^fO, 



i = 0,l,2, , n— 2, 

>7o = 0, 

The theorem involved here may be stated as follows : 

Given f (x, y, v), g (x, y, v) real continuous functions of the three real 
variables (x, y, v ) in the region R: \x — x \<a, \y — y \ < b, \v — v \<b, 
satisfying the inequalities 

\f{x,y',v') — f{x,y,v) | <h 1 \y'—y\+lc 2 \v'—v\, 
1 9 ( x > V'> v ') —9 (x,y,v)\<W y'—y | +&, | v'—v j , 
k lf h 2 being constants and (x, y', v'), (x,y,v) being points in R; then there 
exist solutions y(x), v{x) of the differential equations 

d £ = f{x,y,v), f x = g{x,y,v), 

satisfying the initial conditions y(x )—y , v(x )=v . These functions are 
approached by the solutions of the difference equation system 

Vi+i— yi—i x i+i— %i)f(%i, Vi, v<), 

v i+1 —v i = (x i+1 —Xi)g(x { , y u v { ), i = l, 2, , n. 
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in which x n+% ^=x is a point of the interval \x — x \ < Jc in which the solutions 
of the differential equations are continuous, x 1} x % , . . . . , x n being points of the 
interval (x x) ; the approach mentioned occurs when the number of divisions 
of the interval (x x) is increased indefinitely and the length of the longest 
interval (cCiX i+1 ) is decreased indefinitely.* 

Furthermore, the solutions of the difference system approach y(x), v(x) 
uniformly in any interval which, with its end points, lies within the largest 
interval in which the solutions y(x), v(x) of the differential equations are 
continuous and /, g satisfy the conditions originally imposed, f 

An examination of the proof of this theorem shows that it holds with no 
essential changes if in the difference equations we replace f(x t , y if v { ), 
ff( x t, Vn ih) by functions which approach these values uniformly for all values 
of x' in the interval of convergence of the solution y{x), v(x) when i increases 
with n in such a way that Lt x t =x'. It is of course essential that the sub- 

scripts in the difference equation system be such that y f , v t can be determined 
in succession when y , v are given. 

It is this generalization of the theorem which is of use in the present 
problem. 

The differential equation ^(u) -\-%u=0 is equivalent to the system 
du 



10) 



dv (P—Q'—Mu—R'v , x 

Tx = - ^R = 9{*,y,v), 



and the difference equation system 9) is equivalent to the system 

yii+i— m = {Xi+i—Xi) v t , 

[P [x { ] —Q> [x t ] -M Vi -R' [x t ] v t 



11) 



v i+1 — Vi=(x i+1 — x t )- 



R[x 4 ] 

which evidently corresponds to the above system of differential equations in 
applying the Cauchy— Lipschitz theorem. The functions /, g satisfy a Lipschitz 
condition since R^pO and /, g have continuous first partial derivatives with 
regard to x, u, v in the neighborhood of the solutions u(x), v(x) satisfying the 
initial conditions u(x ) =0, v(x ) — y. 

But instead of the difference system 9) we should like to be able to sub- 
stitute system B, Part II. More precisely, we wish to show that for 57 sufficiently 
large, solutions of B exist uniformly close to corresponding solutions of system 
8) for certain values of X, 

* Pioard, "Traits d'analyse," t. 2, p. 322. 

f Picard, "Traits d'analyse," t. 2, p. 332, and Lipschitz, "Lehrbuch der Analysis," Bd. 2, p. 504. 
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System B is equivalent to 



12) 



Vi+i—Vi— («*+i— a><) Vt , 



v 



i+X 



-v. 



(*<+!-«<) 



(fi + l- % <?<+1 -^.+ ^»'«)-f S<hr7{<fl 



+"Qi+2— Qi+l)^ 



#i+2+A& +2 

i = 0, 1, 2, . . . ., w — 2. 

From the extension mentioned of the Cauchy-Lipschitz theorem we see that 
we can replace system 11) by system 12) if we can show that, as n increases, 
P i+l , (Q i+2 —Qi+i)/&, (R i+2 — R i+1 )/A approach uniformly to P[x { ], Q'OJ, 
B'[««]. 

That P i+1 approaches uniformly to P[x i+1 ] has already been shown. Then 
P i+1 approaches uniformly to P[x { ] because of the uniform continuity of P[x] 
in the interval aj <#<#,»• 

To show that (Q i+Z — Qi+i)/A approaches uniformly to Q'{%i~\, we note that 
Qt+»— Q i+ i= (»«+,— *«+i)G,(*+l + 0i) + (Y i+z —Y i+1 )Q y (i+l-\-d z ) 

in which (*+l + 0) means that the arguments are between (i+2) and (i+1). 
Dividing through by x i+2 — x i+1 , «. e. by A, and using the information regarding 
limits including equation 6), we see that Lt of the right-hand side exists. 

n=oo, t+2=oo 
Qi+Z Qi+1 



Lt 

«=oo, *'+2=8 



%. e. 



QA^+il -\-y'(os i+1 )Q y [x i+1 ] +y"{x i+1 )Q p {x i+1 \ ; 



Lt (Q i+ »-Q i+ i)/A = Q'[x i+1 \. 

n=<xi, i+2=oo 



Furthermore, the limits on the right are approached uniformly, and from this 
fact, together with the uniform continuity of y{x), y' (x) and Q'(x) in the 
interval x <x<x n , it follows that (Q i+Z — Qi+i)/A approaches uniformly to 
Q'lXi], i = l,2, . . . ., n — 2, when n increases indefinitely. Exactly similar reas- 
oning gives similar results for (R i+2 — R i+1 )/A. 

Thus the solution of system 8) is approached by the solution of the 
difference equation system 



13) 



R i+ i s — - + -£- • -£ + [~^- -Pi+^jm 



+ Q 



i+i 






0, 



•=1,2, 



n- 



. Vo = 0, vj 1 = yA. 
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Furthermore, since the coefficients P, Q', R' in the differential system 10) are 
of class C and F pp zp0 in a region of the xy plane including the locus y=y(x), 
%Q<%<% n and for all values of the third argument, a solution of system 11) or 
12) approaches a corresponding solution of system 10) uniformly in the 
interval sc <#<#„, when w increases indefinitely; and similarly the solutions 
of system 13) approach uniformly to the solutions of system 8) as» increases 
indefinitely. 

Then for ^=X. — h and n sufficiently large the solutions of 13) are always 
positive, i>0 and Yi n (% — §)>0, while for A=^, -)-5 and n sufficiently large, 
vt n (% +$)<0. Now choose n=N sufficiently large that y] n (X —S)>0, vj n (X o +S)<0 
and hold n fast. Then >7„(X) is a continuous function of /I, and for some W 
in the interval 1 — $<%<% -\-8, and therefore negative, we have vi n (X)=0. 
Thus, if there is a solution of ^(m) = 0> vanishing at x=x and again at x=x', 
where x <x' <x n , and not vanishing identically, then there exists a W less than 
zero such that in the corresponding solution (*7 >7i • • • • »7»-i>7»)> not identically 
zero, the variable y; n is zero. But this contradicts condition III', and therefore 
a solution u(x) of ^(w)=0 such that u(x ) = and u(x)=£0 can not vanish 
in the interval x < x < x n . 

The results of this work can be summarized in the following 

Theorem I. // the following conditions hold, 

1) {Y x Y t . . . ■ Y n _ 1 ) such that {x x Y x ),. . . .,(*„_il'„_i) lie entirely in the 
interior of region R give a relative minimum for ^>(«/i2/ 2 • ■ • • y n -\) 
in R; 

2) As n increases indefinitely {Y 1 Y % . . . . Y n _ x ) tend uniformly toward 
the corresponding ordinates of a curve y(x) of class C" lying entirely 
in the interior of R and going through the fixed points (x , y ), 

3) When j is allowed to increase with n in such a way that Xj approaches 
a definite value x', Lt (Y s — Y f _{)/A exists and is equal to 

71=00, 7*= CO 

dy(x)/dx\ x=x ,; 

4) F pp [x, y(x), y'{x)]4=0, % <x<x n (y'(x) = dy(x)/dx) : 
then it follows that 

I. y{x) satisfies Euler's differential equation, 

F y [x,y(x), y'(x)]—^F p [x,y{x), y'(x)]=0. 

II. R>0 along* y(x) ; i.e., F pp [x, y{%), y' (<*>)] > 0, x <x<x n . 

* Notice that under hypotheses 1, 2, 3 we have proved that F PP [a>, y(oo), !/'(#)] >0. 
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III. If u(x) is a solution of Jacobi's equation, 

^L{Bu>)~{P-Q')u=Q, 

vanishing for x=x but not vanishing identically, then u(x) does not vanish 
in the interval x <x<x n . 

Thus we see that y (x) satisfies the first three necessary conditions * 
usually given for a function giving a relative minimum of $% n a F(x, y, y')dx. 

It is well to notice that we can make this theorem somewhat stronger by 
replacing the first hypothesis by the part of it which we have really used; 
namely, instead of assuming that (Y x . . . . Y n _ x ) furnishes a relative minimum 
for <p in B, we need merely assume the three necessary conditions for a relative 
minimum in the algebraic problem. 

We note also that our conclusions are not sufficient conditions for y(x) to 
give even a weak minimum of f%«F(x, y, y')dx. Furthermore, if we assume 
sufficient conditions for a proper minimum in the algebraic problem, and 
nothing more, we should not be surprised at being unable to show that the 
conjugate point f x' is beyond x n , for the form A may be positive for any par- 
ticular n and still approach zero as n increases indefinitely. However, if we 
assume sufficient conditions for a proper minimum in the algebraic problem and 
assume furthermore that as n increases indefinitely the least value X of A 
approaches a positive limit, not zero, then our method shows that the conjugate 
point x' is beyond x n if existent and our conclusions are sufficient conditions 
that y (x) give a weak minimum for f^F(x, y, y')dx. 

The above discussion leads to the consideration of the following 

Theorem II. If we assume, 

n 

1) (¥ 1 . ', . . . Y n _ 1 ) give the absolute minimum for ^F[x i} y if (y { — 
2/*-i)M]A, (Xi, y t ) being in region B; 

2) As n increases indefinitely, (Y 1 . . . . Y n _ x ) approach uniformly the 
values of the corresponding ordinates of a curve y=y(%) lying in 
the interior of B, of class C '; 

3) As j increases with n in such a way that Lt x t =x' , 

Lt (Yi-Y^/A 

7l=oo, ;=oo 

exists and is equal to y'(x'), and this limit is approached uniformly 
in the interval x <x'<x n : 

* See Bolza, "Lectures," pp. 101, 102. 

f For definition of conjugate point, see Bolza, "Lectures," p. 60. 
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then it follows that 

I. Lt 2 F[x tt T tt (Y i -Y i _,)/A]A=^F{x, y(x), y'(x)]dx. 

II. // £7iere is a curve of class C lying entirely within R and giving 
the absolate minimum for f^F(x,y,y')dx, then y(x) is such a curve. 

Proof: From hypotheses 2), 3) and the uniform continuity of F(x, y, p) 
for (x, y) in R and p varying between any two fixed limits, it follows that 
as n increases indefinitely F\x t , Y i} (Y { — I r i _ 1 )/A] approaches uniformly to 
F[Xi, yiXi), y' (x { ) ], i=l, 2, .-..., n. From this in turn follows conclusion I. 

Let y x . . . .y n - X be the ordinates of any other curve y—y(x) corresponding 
to the abscissas x x , . . . . , x n _ x , y(x) being of class C and y=y(x) being interior 
to region R. We wish to show that 

f%F(x, y(x),y'(x))dx>f%F(x, y, y')dx. 

n 

For w sufficiently large 2 F\_x i , y if (y t — y i _ 1 )/A]A is as near as we please to 
f%F(x, y, y')dx. But from the first hypothesis 

£ F (x t , y< , ^~^- x ) A > £ F (x t , T t , ¥i ~/ i ~ 1 ) A. 

Then using conclusion I, we read from this inequality conclusion II. 

III. From the Transcendental Problem to the Algebraic Problem. 

Although sufficient conditions for a solution of the algebraic problem do 
not give a solution of the limiting transcendental problem, it is interesting to 
note that if, conversely, y(x) satisfies, sufficient conditions for a weak minimum 
of $%*F(x, y, y')dx, then for n sufficiently large we can find a solution of the 
algebraic problem. To prove this we consider first the 

Theorem : If the following conditions hold, 

1) y=y(x) is a solution of Euler's equation, of class G" in the interval 
x ^x<x n , going through the points (x ,y ), (x n , y n ) and lying 
entirely in the interior of the region* R; 

2) The conjugate point to x is beyond x n or non-existent; 

3) F pp is positive along y(x) ; 

4) y(x) is of class C in an interval x — 5 < a; < rr„ + $, where 8 is some 
positive constant: 

then, given e>0, N exists, independent of i, such that whenever n>N a solution 
of the system 

* The notation here is the same as at the beginning of Part II. 
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F p (i+1)-F p (i) 
A 



-0, 



i=l,2,. 



n- 



1, 



exists, satisfying the relations 



a) 



Vi /'- 1 -y' (x t ) 



<e, 



:1, 2, , n—1, 



A 

») yn=y(x n )- 

We shall show first that for n sufficiently large the system 

Fyii) -EAi±±y=^> 

yo=y(®o), 



r. 



= o, 



1,2, 



w — 1, 



yi—y 



— y' («o) 



has a unique solution satisfying conditions a). 

Let x, y, p be taken as the rectangular coordinates of a point in space. 
Since F pp is positive along y(x), we can construct a "tube" S including as 
interior points all the points [x, y (x) , y' (x) ] and such that F pp (x, y, p) is 
positive for all points of this tube.* Let 

1) v=F p (x, y+p-A, p) 

define a transformation of the points of the four-space (x,y, p ; A) such that 
(x, y, p) is a point in 8 and A is a small positive quantity, into points of 
(x, y, v; A) space, and let L denote the image of y{x), i.e., of [x,y(x), 
y'(x) ; 0]. We shall denote the locus |>, y{x), y' (x) ; 0] by k. The function 
F p (x, y-\-p-A, p) is of class C in a region of (x, y, p; A) space including as 
interior points the points of the locus k. Differentiating this function with 
respect to p, we have 

Fpp(x, y+V • ^ p) +A-F py (x, y+p -A,p), 
which is positive along the closed interior locus k. All the hypotheses of the 
extended theorem of inversion f are fulfilled, so that 

a) A positive quantity p can be found so small that to each point 
(x, y, v; A) in (L, p), the image of the p-neighborhood of k, k p , there corre- 
sponds one and but one point (x, y, p; A) of k p . 

b) The inverse function $, 

2) P=<p{oc, y, v; A), 
is of class C in (L, p). 



* See Bolza, " Vorlesungen fiber Variationsrechnung," p. 157, Lemma II. 
t See Boljsa, "Vorlesungen," p. 160, b; p. 163, c. 
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c) Corresponding to p we can find a positive quantity <r such that to each 
point of the a-neighborhood of L, L„ there corresponds, through the trans- 
formation 2 ) , one and only one point in h p . 

If there is a solution of I' satisfying conditions a) as desired, choose A 
so small that [x t y, (y t — y^x) /A; A] lies in fe p , i = l,2, ...., n — 1. Let 
F p (x iy y«_i + Ap<, p^ = v { _ lf where p«=(y«— y*_i)/A. Then from 2) we have 
Pi=<p(x t , yi_i, Vi_ x ; A), and the solution of system I' lying in h p becomes a 
solution of the following system lying in (L, p) : 

Vi— Vi_i = (*<+i— a?i) F, S a>« , 2/i_! + A • <J> (a;, , y,_ x , «<_, ; A) , 
4>(^,«/»-i, v f _ l5 A) jsA-/j(aj ?/,_!, v^x; A), 

y =y(«o); (yi— y <>)/&= y'(x )- 



II. 



Conversely, if II has a solution in L^ , then I' has a solution in h„ . Thus the 
system II is equivalent to system I'. 
Euler's differential equation, 

F y (x,y,y')--^F p (x,y,y f ) = 0, 

is equivalent to the system 

dv 



III. 



dx =F y [x,y,®(x,y,v)]^g(x,y,v), 
■g| =<&(*, 2/, v)*EBf(x, y, v), 



in which <t>(a, ;?/, v) = $(x, y, v ; 0). Then our hypothesis is that !>(#),:?/(#)] 
is a set of solutions of III, y (x) being the original extremal and v (x) = 
F p [x,y(x), y'(x)]. The functions g (x, y, v), f (x, y, v ) satisfy a Lipschitz 
condition in the sub-region of (L, p) in which A=0, since they are of class C 
in this region. Then the solution [v(x), y(x)] of III is approached by a 
solution of the difference equation system, 

, Vi— v t _! = (x { — a^-i) F y [>,_, , y t _ x , O (^_ x , y t _ x , «,_!> ] , 



Vi— tft-i= ( x i~ **-i)*(**-n 2/t-i, «<-i), *=1, 2, , n— 1. 

Comparing this system with II, we note that since <p{x, y, v; A) is of class C" 
in (L, p) and F y is continuous in its three arguments in the region of interest, 
the functions on the right in II approach uniformly to the corresponding 
functions on the right in IV as n increases indefinitely, t=l, 2, ....,n — 1. 
Then the generalization of the Cauchy— Lipschitz theorem, specified in Part II, 
informs us that for n sufficiently large system II has a solution as near as we 
48 
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please to [v (x), y(x)] ; i. e., given e>0, N can be found such that whenever 
n>N, 

I Vi—Vi^i) | < e, I v t —v (Xi) | < e. 

Furthermore, N can be chosen independent of i. Here we use the fact 
that the Cauchy-Lipschitz method gives, for n sufficiently large, solutions of a 
difference equation system uniformly close to the solution of the differential 
equation system in any interval which, together with its end points, lies in an 
interval in which the solutions of the differential system are continuous. We 
have assumed y{x) to be of class C" for x — 5<#<;c„+^ 5>0. That v(x) 
is continuous in such an interval follows from this restriction on y(x) and the 
continuity of F p for (x, y) in R and for all values of p. 

We have, also, 

yi—yi-x _ y,^ _ ^ (g. o y._ x f V{ _ i; A)—(p(Xi, y(Xt), v(x t ); 0). 
Since $ is continuous in (L, p), N exists independent of i such that 



y±^=l-y'{ Xi ) 



<6, i = l,2,. 



whenever n>N. It is evident that the solutions of systems II and I' are unique. 
Having established the existence of a unique solution of I' satisfying con- 
ditions a), we are in a position to use the idea of a set of extremals through a 
point to prove the theorem stated above. Since we assume that the conjugate 
point of x is beyond x n or non-existent, y{x) is one of a set of extremals* 
y(x, y) through (x , y ), y being the slope of the curve y=y(x, y) for x = x , 
which do not intersect in the interval x <x<x n . Let us write y = y'(^o)- 
Then there exists a positive number d x such that when | y — y 1 < d x the corre- 
sponding extremals do not intersect in the interval x <x<x n . Also a positive 
number d 2 exists such that when \y — yo\<d t , F Pp is positive along y(x, y). 
Then the hypotheses for the application of the Cauchy-Lipschitz theorem are 
satisfied in the neighborhood of any one of the curves y (x, y'), where | y<>~y' \ <d, 
d being the smaller of a\ and d 2 . Furthermore, an examination of the Cauchy- 
Lipschitz proof shows that the theorem applies when the initial value of v, v 
is made a function of the parameter y, \y — y\ < d, and therefore we have that 
for a given e > 0, N t can be found independent of i and y such that, whenever 
n>N iy 

y{(y)—y(«>i,r)\ <<?, 

yi{y) ~l i -' {r) -y'{*i,r)\<e, *=i, 2, ....,»; | r _ yo |<<*, 

in which the set y x {y), . . . ., y n {y) is the solution of I' for y'(x ) =y. 

* See Bolza, "Lectures," p. 61. 
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Our conclusion is that for n sufficiently large we are able to get the solu- 
tion of the difference system I' uniformly close to the corresponding extremal 
throughout the interval % <a;<% n and for \y — y \ <d, y being the initial slope 
of the extremal in question. We are now in a position to prove the original 
theorem ; *. e., to find a solution of the system I satisfying conditions a) and b) . 

As a function of y alone, y(x, y) and y'(x, y) are continuous for 
\y — y \ <d. Then a positive number d B exists such that 

\y (*,y)— y foyo) I <*/2, 

\y'(x, y)— y'(x,y )\ < e/2, 

whenever \y—y \<d s . Let y lt y 2 be such that y!<y <y 2 , y 2 — yo<^ 3 , y — 
yi<d 3 . As a function of y alone, y(x, y) is an increasing function, since the 
set of extremals through (x , y ) do not intersect in the interval in question, 
and there exists a positive number ft < e/2 such that 

y(®«, yd—y^n, y )>h, 
y(v n , yo)—y(Xn, yi)>h- 

Now choose N 2 so large that for n>N 2 

\yikr)—y(^i,r)\<K 
ydy) —yi-i(r) 

A 

Choose n—n'>N 2 and hold it fast. Then the inequalities y„>(yi) <y n '(yo) < 
yn'(yz) hold. Also y n >{y) is a continuous function of y, \y — y \ <d % . Then 
for some value of y, y=y' such that yi<y'<y 2 ) y n '(y') =y( x n> yo)- Further- 
more, tyi(y') — y(x if y') | <e/2 and \y(x t , y')—y(Xi, y ) | < e/2. 

••• \yAy')— y(®i, yo) I <*, 



y'(<*>i, y) 



<h, 



yi<y<y 2 . 



and similarly, 



yi(y')—yi~i(y') 

A 



<s. 



y'(®i, yo) 

, y„'(y') we have a solution of system I satisfying con- 



Thus, in yi(y'), . 
ditions a) and b). 

To carry out the plan of Part III we wish now to show that for n suffi- 
ciently large a minimum of the algebraic problem exists. We shall show that 
(y lf . . . ., y n _ x ) exist satisfying the sufficient conditions I', II' III" of Part II 
for the desired minimum. We already have the existence for n sufficiently 
large of a solution (y lf . . . . , y n - X ) of 

FJi+l)-FJi) _, 



F y {i) 



0, 



i=l,2, . . . ., n — 1, 
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going through. (x , y ), (x n , y n ) and as near to the extremal y(x) as we please. 
To show that, for n sufficiently large, condition 

II. A«P«+2Ag,+fl, +1 +fl,>0, i=l, 2, ...., n-1, 

is satisfied, we note that we can make R it R i+1 as near to F^x^ y{x t ), y' (x t )] 
as we please and therefore positive, and that P i} Q { remain finite. It is sufficient 
now to show that the difference equation system B, Part II, has no solution 
(>7o > *7i > • • • •> V n ) n °t vanishing identically, such that yj =yi n =0 for X,<0. (It 
should be noticed that the equality sign in X<0 is especially important.) 

We have in Part II that for n sufficiently large a solution of the difference 
system B such that yi =§, (yi 1} . . . ., v n -i)^^ is uniformly close to a corre- 
sponding solution of the system 

d 



3) 



^(I2>/) + (<?'-P+;i)>7=0, 



y 1 (x ) = 0; vj(x)^0. 

But by hypothesis the conjugate point to x is beyond x n or non-existent. This 
means that if solutions of 3) for ^,=0 vanish for any x greater than x , the 
next root, x', is greater than x n ; then a theorem* of Sturm's tells us that for 
X < the next root would be still further to the right. It is therefore evident 
that for n sufficiently large no solution of the difference system exists for A<0 
such that ^ =>7„=0 and {yj x , . . . ., >7„_ x )^0. Thus for n sufficiently large we 
have the existence of {y x , . . . . , y n -x) satisfying the sufficient conditions stated 
in Part II for a relative minimum of the algebraic sum. Furthermore, it is 
well to note that for n sufficiently large the ordinates (y t , . . . ., 2/„_i) are as 
close as we please to the extremal giving the assumed minimum of the trans- 
cendental problem. 

IV. The Tkansobndental Pboblbm Following Foemal Analogies 

TO THE AlGEBBAIC PbOBLEM. 

There is evidently a very close correspondence between the algebraic and 
the transcendental problems which we have considered. The purpose of this 
part of the paper is to make this correspondence still more striking by bringing 
out the formal analogy involved. 

Let F(x, y, p) be a function of class C" in a certain simply connected 
region R of the xy plane including as interior points the points (0,0) and (1,0), 
and for all values of p. Let y (x) be a function of class C" in the interval 
0<.-r<l, such that the curve y—y(x) lies entirely in the interior of R and 
goes through the points (0, 0) and (1, 0). Let 

* See Theorem b) , Part II. 



{ 
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i(y) = Po(a?)^+Pi(«)j|+p 8 (*) -y 

be a self-adjoint differential expression such that the system 

L(y)=0, 
y(0)=y(l)=0 

is incompatible, and let G (x, £) be the corresponding Green's function.* Suppose 
we note that L[y{x)~] = r(x). Let yj{x) be continuous in the interval 0<r»<l, 
and let y(x) be the function determined by 

L(y) = r(x) + eY{(x), 

y{0)=y{l) = Q, 

where e is a constant chosen small enough that y(x) lies in the interior of R. 
Making use of the Green's function, we may write 

y(*)=Shr(e)G{x,Z)dS; y'(x)=ftr(Z)G 1 (x,£)dl; > 

where G x {x, £) = ^G(x,£); 

y{x)=Sllr{Z)+BY,{Z)]G{x,Z)dZ; y'{x) = $Hr{Z)+eYi{Z)}G x {x,Z)dZ; 
a) y{x)-y{x)=e$lY,{Z)G(x,!r)dZ; y' (x) -y'(po) =e$\n{Z)G 1 {x,Z)dZ. 
Using Taylor's series with a remainder, 
F{x,y,y') = F{x,y,y') + (y-y)F y +(y'-y')F p 

JL< — 

2! 

3! 

the arguments of all partial derivatives except in the last group being (x,y,y'), 
those of the last group being of the form [x, y + 6 x (y — y), y'-\-O z (y' — y')], 
d lt 2 being proper fractions. Substituting the values of y — y, y' — y' from a), 
we have 
F{x,y,V) = F{x,y,y')+e\FJh n {Z)G{x,Z)dZ+FJU{Z)G x (x,Z)dZ\ 

+ (- l \Fy y (nv^)G(x^)dO 2 +2F vp Sly 1 ^)G(x,^d^S 1 yia)Gi(^)d^ 

+* w (JSi(«fl ! i(»,*)«i*) , | +«•(*) 

= F{x, y ,y')+efa{Z)[G(x J Z)F y +G 1 {x,Z)F p \dZ 

+ Y^Poyi{^)n{^){G{x,^)G{x,^)F yy +2G{x^ l )G 1 {x,^)F yv 

+ G 1 (x,Z 1 )G 1 (x,&)F pp }dZ 1 dZ 2 +6*(e). 

* For the definition of the Green's function see an article by B6cher, Annals of Mathematics, 
Vol. XIII, p. 71. 



+ ^\(y-yyF yy +2(y-y)(y'-y>)F yp + (y'-y>yF pp \ 

+ ^ri (y-y) s F yyy + . . . . + (y'-y') s F ppp \, 
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Suppose y(x) makes flF(x, y, y')dx as small as or smaller than does any- 
other curve of class G" in a certain neighborhood of y = y (x) and going 
through the points (0,0) and (1,0). Then for e sufficiently small 

Fo[F(x,y,y')-F(x,y,y')lda> = 0. 
Employing reasoning which was useful in Parts I, II, we obtain conditions 
exactly analogous to conditions \here derived, 

H". Sl\Slftn(Zx)n(&) [G{x,^)G{x,^)F yy +2G{x,^)G 1 {x,^)F w 

+ G 1 {x,Z 1 )G x {x,&)F pp \d&d&\dx>0. 

Up to this point we have made no essential use of the fact that the third 
argument in F is the derivatives of y (x), or in any way related to y (x). We 
should have come to conditions analogous to I" and II" if the third argument 
in F had been any other function, of x and had been varied similarly to the 
way in which y(x) was varied. So far the important fact has been the possi- 
bility of expansion by means of Taylor's series with a remainder, i. e., an 
expansion in ascending powers of a quantity e which can be made as small as 
we please. From the very definition of maximum and minimum it is evident 
that such an expansion is fundamental in a general analytic investigation of 
such problems, whether the variable be a point or a curve. 

Since G(x,%) is continuous in the square 0<a?<l, 
0<£<1> ari d G x (x,^) is continuous in this square except 
along the diagonal x=%, where it has a unit discontinuity, 
the order of integration may be reversed in I" and II", 
giving 

I'. fad) [fo\G{x,Z)F y +G x (x,Z)F p \dxW=Q, 

ir. tf/u (&)»?<&) [J3 W*,*i)0(*»*.)-p 

+ 2G(x,£ 1 )G 1 (x,&)Q + Gi(oc,Si)G 1 {x,Z z )R\dx]dZ 1 dZ 2 

where 4> (£i , £ 2 ) represents the integral in brackets [ ] . 

Since I' is an expression for the first variation of flF (x,y,y')dx, we 
should expect to be able to show from it that y(x) is a solution of Euler's 
equation. To do this we proceed as follows : 

ftGi(*,£)Fpdx = ftG 1 (x,Z)F,dx+ftG 1 (x,£)F v da>. 
Since in each sub-interval G x {x,%) is continuous, we may integrate by parts 
thus : 

SWi(*,Z)F p dx=G(x,Z)F p \t-ftG{x,iZ)^(F p )dx, 

ftG 1 (x,Z)F p dx=G(x,!;)F p \\-S\G(x,j;)-^(F p )dx. 




Bobbins: A Method in the Calculus of Variations. 391 



Since G(x,£) is continuous and 6r(0,£) = G(1,%) =0, 

S 1 oGi(x,£)F P dx = -f G(x,S)^(F p )dx. 
Then I' may be written, 

Since 1) holds for yi(%) any continuous function, it holds if >?(£) is of class C 
in the interval 0<£<1 and >?(0) = vi(l) = 0. The coefficient of *}{%) in the 

integrand involved in 1), flG(x,%) \F y — -=-F p \dx, is a continuous function 

of §. Then, by the fundamental lemma* of the Calculus of Variations, 

A 
dx~ 



r G(x,£)\F y -dF p \dx^0. 



From this it follows that 

I. F y {x,y{x),y'{x)^— -^F p [x,y(x),y'(x))^0. 

For, suppose F y ~-j- F p = q>( x) ^0. Since G(x,g) is symmetric, we have 
$\G{x, %)q>(x)dx = 0; but this would be the solution of the system 

u(0)=u(l)=0, 

which can not vanish identically since $(o;)=£0. The contradiction proves 
the truth of I. 

What we have really proved is that of all curves of class G" through the 
points (0,0), (1,0), those giving a minimum must satisfy Euler's equation. 
We can not expect more of this method since we consider no more general 
curves. However, we shall assume in the discussion of the second variation 
that F pp =f=0 along y{x) ; and it has been shownf that in this case an extremal 
of class G' is of class G" . 

It is interesting to notice that the first variation, I', is a simple integral 
with respect to £; that the second variation, /J/o>7(§i)>7 (&!)<?>(£i> ^d^d^ 
can be considered as a double integral ; and it is evident that the k-th. variation 
would be a Jc-iold integral. Here we have a complete formal analogy to the 
simple and multiple sums in the corresponding variations in the algebraic 
problem. 

In examining the second variation, II', we shall prove first that 4>(£,£)>0, 
0<£<1. This condition is analogous to the condition that $«>0 in the algebraic 
problem, $ and the method of proof is very similar. 

* See Bolza, "Lectures," p. 20. f Bolza, "Lectures," pp. 25, 26. $ See Part I, condition 2). 
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Suppose that for £=%', $(£, £) = — c, where c is positive, x' being a point 
in the interval 0<a;<l. The function 4>(£i, £ 2 ) is continuous in the square 
O^i^^l, 0<£ 2 <1, as can be seen by breaking up the integral representing 
4>(£i> £2) i n to integrals over intervals in which the integrand is continuous. 
Then we can find a positive quantity h such that 4>(§i, £ 2 ) <0 in the square, 
x' — ^<£i<#' + ^> x' — fc<£ 2 <#' + &. Now choose yi{x) positive in the interval 
x' — h<x<x'-\-h, zero outside of this interval, and continuous in the interval 
<o? < 1. Then we have the inequality 

which is contrary to condition II'. 

■•.*<*,*)>o, 0<£<1. 

From this condition we shall show that if R(x)=f z in the interval 0<o;<l, 
then R(x) is positive in this interval. For this purpose we need the following 
properties of G(x, £), G 1 (x,^) : 

1) G(x,%) as a function of x satisfies the boundary conditions u(0) = 
«(1)=0; i.e., G(0,£) = <?(1,£)=0. 

2) G(x,%) is symmetric. .'. G(x,0) = G{x,l)=0. 

3) G 1 (x,^) x=it — G 1 (x,%) x = i _=l. (x=$- + is read, "as x approaches £ from 
the right.") 

These facts enable us to show by direct computation that <jb(0,0) =0 and 
4»'(0,0)=B(0), where 4.' (£, £) = dq> (£, Q/&Z,. We see that $'(0,0) is not 
negative, since if it were, $(£, £) would be negative for very small values of £. 
.*. 22(0) >0, since we assume R(%)zfzQ in the interval 0<£<1; then the con- 
tinuity of R(%) gives us the desired result, 

II. R(x)>0, 

in the interval < x < 1. 

In the further consideration of the second variation we shall show that as 
the kernel* of an integral, 4>(^> £ 2 ) has only positive characteristic numbers. 
Suppose a,=3. is a characteristic number and that 

»7 (») = *<> .ft $ (a;, £)>7 (£)«*£• 
Multiplying both sides by vi(x) and integrating, 

SW{x)dx = \§l§l<p{x,^)ri(^)Yi{x)d^dx. 

Since the left-hand side is positive and the double integral is not negative, 
and therefore actually positive in this case, X must be positive. 

* The number X is a characteristic number of a kernel F (x, y) if a continuous function u (x) exists 
such that u(tx>)=\f <1 1 F(x, |)m(?)<2£. See also any book on integral equations. 
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To get Jacobi's condition we proceed as follows: Let y{x) — y{x) = eu(x), 
where u(x) is of class C" and «(0) = «(1) = 0. It is well known* that for 
such a variation on y (cc) the second variation of the integral can be written 

2) J(u) 33 p o wi>{u)dx, 

wheref ^(u) = (P— Q')u— d(Bu')/dx, i. e., 

J(u) = Sl[(P—Q')u*--u-^ ; (Ru')]dx>0. 

It is now convenient to consider only such variations u(x) that §\u 2 (x) dx — \. 
This restriction is exactly analogous to those made in Parts I, II in dealing 

with the second variation in the algebraic problems, the restrictions 2 v\ — 1, 



4=1 



2 rft A=l, and it should be noted that in no case is this an essential restriction. 



i=l 



Then, if s is a positive constant, J(u) +e is positive: 

j( M )+ e== j-i[(P_ Q' + e )^_ u~{Ru')]dx>0. 

Introducing the notation 

L(u)^^(Ru') + (~Q'-P-e)u, 

we have 

— fluL(u) dx > 0. 

Hence, we conclude that the system 

L(u)=0, 



3) I u(0) = u(l)=0 

is incompatible, for if it were compatible, u(x) being a solution, a positive 
constant c could be found such that $lu?dx = l and §luL(u)dx — 0, where 
u — cu. Let G(x,%) be the Green's function of system 3), and let the system 

f L(u) = yi(x), 
{ u (0) = u{l) = 0, 

ri{x) being a continuous function, 0<#<1, define a variation of the curve 
y(x) ; (i. e., a function u such that y — y=eu). Then 

and 

fluL(u)dx = /J /Jj7(£) »?(«) <?(*>,£) df d» < 0. 



* See Bolza, "Lectures," p. 52. 

f \f/(u) as here used is a standard notation; it should be noted that it is the negative of the f(u) 
used in Part II. 

49 
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The characteristic numbers of the kernel G(x,%) are all negative. For suppose 
that one, X > were positive, and that 

n{x)=^$lG(x,Z)Yi(Z)d£. 
Multiplying both sides by yi{x) and integrating, we have 

flyi*(x)dx=l Stfty ! (£) v (x)G(x,£)d£dv. 

The left-hand side is positive, while the right-hand side is negative, and the 
contradiction proves the statement. But the characteristic numbers of the 
kernel G(x,%) are the same as for the system* 

L(u) = Xw, 
w (0) = w(l) = 0. 

Then the system 

,, [ L(u) + lu=Q, 

' l«(0) = «(l)=0 

has a solution u(x)^0 only for positive values of X, no matter bow small 
a positive constant e be chosen. Then the system 



{ 



5) 



A c (Ru') + (Q'-P)u + *u = 0, 



!*(0)=«*(1)=0 

has no negative characteristic numbers, for if Jl < were a characteristic 
number and we choose e< — 7i , then 4) would have a negative characteristic 
number, namely, X, +«. 

There is one and but one valuef of X for which system 5) has a solution 
not vanishing in the interval 0<o;<l, and we have shown that this value is 
positive. Then there is no non-identically vanishing solution of the system 

-4(u),= -^(Ru') + (Q' + P)u=0, 

w(0)=0; u(x')=0, 0<cc'<l, 

for if there were, system 5) would have a solution vanishing nowhere in the 
interval < x < 1 for "K negative. % 

* See, for instance, Hilbert, "Grundziige einer allgemeinen Theorie der linearen Integralgleiehungen," 
Oottingische Nachrichten, 1904, p. 225. 
t See Theorem a), Part II. 
f See Theorem b), Part II. 



